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11 Summary. 

12 1. Understanding space usage and resource selection is a primary focus of many studies 

13 of animal populations. Usually, such studies are based on location data obtained from 

14 telemetry, and resource selection functions (RSF) are used for inference. Another important 

15 focus of wildlife research is estimation and modeling population size and density. Recently 

16 developed spatial capture-recapture (SCR) models accomplish this objective using individual 

17 encounter history data with auxiliary spatial information on location of capture. SCR models 

18 include encounter probability functions that are intuitively related to RSFs, but to date, no 

19 one has extended SCR models to allow for explicit inference about space usage and resource 

20 selection. 

21 2. In this paper we develop the first statistical framework for jointly modeling space usage. 



22 resource selection, and population density by integrating SCR data, such as from camera 

23 traps, mist-nets, or conventional catch-traps, with resource selection data from telemetered 

24 individuals. We provide a framework for estimation based on marginal likelihood, wherein 

25 we estimate simultaneously the parameters of the SCR and RSF models. 

26 3. Our method leads to increases in precision for estimating population density and 

27 parameters of ordinary SCR models. Importantly, we also find that SCR models alone can 

28 estimate parameters of resource selection functions and, as such, SCR methods can be used 

29 as the sole source for studying space-usage; however, precision will be higher when telemetry 

30 data are available. 

31 4. Finally, we find that SCR models using standard symmetric and stationary encounter 

32 probability models produce biased estimates of density when animal space usage is related to 

33 a landscape covariate. Therefore, it is important that space usage be taken into consideration, 

34 if possible, in studies focused on estimating density using capture-recapture methods. 

35 Key- words. animal movement, animal sampling, encounter probability, hierarchical 

36 modeling, landscape connectivity, marginal likelihood, resource selection, space usage, spa- 

37 tial capture-recapture. 

38 

39 1 Introduction 



40 
41 
42 



Spatial capture-recapture (SCR) models are relatively new methods for inference about 
population density from capture-recapture data using auxihary information about individual 
capture locations (lEffordl 120041 iBorchers and Effordl 120081 iRoyle and Youngl 120081). SCR 



models posit that N individuals are located within a region denoted S. Each individual 
has a home range or activity area within which movement occurs during some well-defined 
time interval, and the center of the animal's activity has Cartesian coordinates Sj for in- 
dividuals i = 1, . . . ,N. The population is sampled using J traps with coordinates for 
j = 1, . . . , J, and encounter probability is expressed as a function of the distance between 
trap location (xj), and individual activity center (sj). While SCR models are a relatively re- 



49 cent innovation, their use is already becoming widespread (Efford et al. , 2009 Gardner et al. 



2010b|a Kery et al. 2010 Gopalaswamy et al. 2012; Foster and Harmsen 2012) because 



they resolve critical problems with ordinary non-spatial capture-recapture methods such as 
ill-defined area sampled and heterogeneity in encounter probability due to the juxtaposition 



53 of individuals with traps (Borchers 2011). Furthermore, unlike traditional capture-recapture 



methods, SCR models allow for inference about the processes determining spatial variation 
in population density. 

Despite the increasing popularity of SCR models, every application of them has been 
based on encounter probability models, such as the bivariate normal distribution, that imply 
symmetric and stationary (invariant to translation) models for home range. While such 
simple models might be necessitated in practice by sparse data, home range size and shape 
are often not well represented by stationary distributions because animals select resources 
that are unevenly distributed in space. Therefore more complex models are needed to relate 
the capture process with the way in which individuals utilize space. 

In this paper, we extend SCR capture probability models to accommodate models of space 
usage or resource selection, by extending them to include one or more explicit landscape 
covariates, which the investigator believes might affect how individual animals use space 
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within their home range (this is what (Johnson, 1980) called third- order selection). We do 
this in a way that is entirely consistent with the manner in which parameters of classical 



resource selection functions (RSF) (Manly et al. , 2002) or utilization distributions (UD) 



69 (Worton, 1989 Fieberg and Kochanny, 2005 Fieberg, 2007) are estimated from animal 



telemetry data. In fact, we argue that SCR models and RSF/UD models estimated from 
telemetry are based on the same basic underlying model of space usage. The important 
distinctions between SCR and RSF studies are that (1) resource selection studies do not 
result in estimates of population density and (2) in SCR studies, encounter of individuals is 
imperfect (i.e., "p < 1") whereas, with RSF data obtained by telemetry, encounter is perfect. 
With respect to the latter point, we can think of the RSF and SCR studies as being exactly 
equivalent either if we have a dense array of trapping devices, or if our telemetry apparatus 
samples time or space imperfectly. A key concept that we must confront in order to unify and 
integrate SCR and RSF data is that we need to formulate both models in terms of a common 
latent variable so that we can make them consistent with respect to some underlying space 
utilization process. As we will explain, this latent variable is the number of times that an 
individual uses a particular region of the landscape over some period of time. 

The modeling framework we develop here simultaneously resolves three important prob- 
lems: (1) it generalizes all existing capture probability models for SCR data to accommodate 
realistic patterns of space usage that result in asymmetric and irregular home ranges; (2) it 
allows estimation of RSF parameters directly from SCR data, i.e., absent telemetry data; and 
(3) it provides the basis for integrating telemetry data directly into SCR models to improve 
estimates of model parameters, including density. Our model greatly expands the applied 
relevance of SCR methods for conservation and management, and for addressing applied and 
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89 theoretical questions related to animal space usage and resource selection. 



2 Spatial Capture-Recapture 



91 A number of distinct observation models have been proposed for spatial capture-recapture 



92 studies (Borchers and Efford 2008; Royle et al. 2009; Efford et al. 2009), including Poisson 



93 multinomial, and binomial observation models. Here we focus on the binomial model in which 

94 we suppose that the J traps are operated for K periods (e.g., nights), and the observations are 

95 individual- and trap-specific counts yij, which are binomial with sample size K and capture 

96 probabilities pij which depend on trap locations Xj and individual activity centers Sj as 

97 described subsequently. The vector of trap-specific counts for individual i, = {yn, . . . , yij) 



98 is its encounter history. A standard encounter probability model (Borchers and Efford, 2008) 

99 is the Gaussian model in which 



log(Pii) = "0 + <^id% (1) 

100 or, equivalently, Pij = Aq exp(— c/f^/ (2cr^)), where dij is the Euclidean distance between points 

101 Si and Xj, dij = ||sj — Xj|| = a/ {sh — xji)'^ + (sj2 — a^j2)^, and ao = log{Xo) and ai = 

102 —1/(2(7^). Alternative detection models are used, but all are functions of Euclidean distance 

103 and so we do not consider them further here. 

104 The primary motivation behind our work is that, in all previous applications of SCR 

105 models, simple encounter probability models based only on Euclidean distance have been 

106 used, with estimation based on standard likelihood or Bayesian methods. These methods 

107 regard the activity center for each individual i, Sj, as latent variables and remove them 



5 



108 from the likelihood either under a model of "uniformity" in which s ~ Unif(5) where S is 

109 a spatial region (the "state-space" of s), or a model in which covariates might affect the 



no spatial distribution of individuals (Borchers and Efford, 2008). The state-space S defines 

111 the potential values for any activity center s, e.g., a polygon defining available habitat or 

112 range of the species under study. 

113 A critical problem with standard SCR models is that the encounter probability model 

114 based on Euclidean distance metric is unaffected by habitat or landscape structure, and 

115 it implies that the space used by individuals is stationary and symmetric, which may be 

116 unreasonable in many applications. For example, if the common detection model based on a 

117 bivariate normal probability distribution function is used, then the implied space usage by all 

118 individuals, no matter their location in space or local habitat conditions, is symmetric with 

119 circular contours of usage intensity. Subsequently we provide an extension of this class of 

120 SCR models that accommodates asymmetric, irregular and spatially heterogeneous models 

121 of space usage. Thus, "where" an individual lives on the landscape, and the state of the 

122 surrounding landscape, will determine the character of its usage of space. In particular, we 

123 suggest encounter probability models that imply irregular, asymmetric and non-stationary 

124 home ranges of individuals and that are sensitive to the local landscape being used by an 

125 individual. 



126 3 Basic Model of Space Usage 

127 We develop the model here in terms of a discrete landscape purely for computational 

128 expediency. This formulation will accommodate the vast majority of actual data sets, as 

129 almost all habitat or landscape structure data comes to us in the form of raster data. Let 
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130 Xi, . . . , x^c* identify the center coordinates of a set of nG pixels that define a landscape. In 

131 SCR studies, a subset of the coordinates x will correspond to trap locations where we might 

132 observe individuals whereas, in telemetry studies, animals are observable (by telemetry fixes) 

133 at potentially all coordinates. 

134 Let z{x.) denote a covariate measured (or defined) for every pixel x. For clarity, we 

135 develop the basic ideas here in terms of a single covariate but, in practice, investigators 

136 typically have more than 1 covariate, which poses no additional problems. We suppose that 

137 a population of individuals wanders around space in some manner related to the covariate 

138 z{x.), and their locations accumulate in pixels by some omnipotent accounting mechanism. 

139 We will define "use of x" to be the event that an individual animal appeared in some pixel 

140 X. This is equivalently stated in the literature in terms of individual having selected x. As a 



141 jDiological matter, use is the outcome of individuals moving around their home range (Hooten 



142 et al. , 2010), i.e., where an individual is at any point in time is the result of some movement 

143 process. However, to understand space usage, it is not necessary to entertain explicit models 

144 of movement, just to observe the outcomes, and so we don't elaborate further on what could 

145 be sensible or useful models of movement. 

146 Suppose that an individual is monitored over some period of time and a fixed number, 

147 say R, of use observations are recorded. Let n(x) be the use frequency of pixel x for that 

148 individual, i.e., the number of times that individual used pixel x during some period of time. 

149 We assume the following probability distribution for the nG x 1 vector of use frequencies: 

n ~ Multinom(i?, tt) 
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150 where tt is the nG x 1 vector of use probabihties with elements (for each pixel): 

exp(a22:(x)) 



TT X 



Ex exp(a22;(x)) 



151 This is the standard RSF model (Manly et al. , 2002) used to model telemetry data. The 



152 parameter a2 is the effect of the landscape covariate 2;(x) on the relative probability of use. 

153 Thus, if is positive, the relative probability of use increases as the value of the covariate 

154 increases. In practice, we don't get to observe {?t,(x)} for all individuals but, instead, only 

155 for a small subset say i = 1,2,..., Ntei, which we capture and install telemetry devices on. 

156 For the telemetered individuals, we assume they behave according to the same RSF model 

157 as the population as a whole, which might be justified if individuals are randomly sampled 

158 from the population. 

159 We extend this model slightly to make it more realistic spatially and also consistent 

160 with standard SCR models. Let s denote the centroid of an individuals home range and let 

161 d{s,x) = ||x — s|| be the distance from the home range center s of some individual to pixel 

162 X, and let n(x, s) denote the use frequency of pixel x for an individual with activity center 

163 s. We modify the space usage model to accommodate that space use will be concentrated 



164 around an individual's home range center (Johnson et al. , 2008 Forester et al. , 2009): 



exp(-aic/(x,s)2 + a22;(x)) 
exp(-ai(i(x, s)^ + a2-2(x)) 

165 where ai = l/(2cr^) describes the rate at which encounter probability declines as a function 

166 of distance, (i(x, s). From ordinary telemetry data, it would be possible to estimate param- 

167 eters ai, a2 and also the activity centers s using standard likelihood methods based on the 
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168 multinomial likelihood (Johnson et al. , 2008). 



169 Note that Eq. |2] resembles standard encounter models used in spatial capture-recapture 

170 but with an additional covariate z{x.). The main difference between this observation model 

171 and the standard SCR model is that the model here includes the normalizing constant 

172 J2x exp(— aic/(x, s)^ + q;2-2(x)), which ensures that the use distribution is a proper probability 

173 density function. Thus we are able to characterize the probability of encounter in terms of 

174 both distance from activity center and space use. Note that, under this model for space 

175 usage or resource selection, if there are no covariates, or if a2 = 0, then the probabilities 

176 7r(x|s) are directly proportional to the SCR model for encounter probability. For example, 

177 setting a2 = 0, then this implies probability of use for pixel x is: 

p(x|s) oc exp(— Q;i(i(x, s)^). 

178 Therefore, for whatever model we choose for p(x, s) in an ordinary SCR model, we can 

179 modify the distance component in the RSF function in Eq. |2] accordingly to be consistent 

180 with that model, by choosing '7r(x|s) according to 

7r(x|s) oc exp{log{p{-x\s)) + a2z{'x)) 

181 As an illustration of space usage patterns under this model, we simulated a covariate 

182 that represents variation in habitat structure (Fig. [T]) such as might correspond to habitat 

183 quality. This was simulated by using a simple kriging interpolator of spatial noise. Space 

184 usage patterns for 8 individuals in this landscape are shown in Fig. |2| simulated with 

185 ai = 1/(2(7^) with a = 2 and the coefficient on z(x) set to 02 = 1. These space usage 
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186 densities - "home ranges" - exhibit clear non-stationarity in response to the structure of the 

187 underlying covariate, and they are distinctly asymmetrical. We note that if were set to 0, 

188 the 8 home ranges shown here would resemble bivariate normal kernels with a = 2. Another 

189 interesting thing to note is that the activity centers are not typically located in the pixel of 

190 highest use or even the centroid of usage. That is, the observed "average" location is not an 

191 unbiased estimator of s under the model in Eq. [2j 

192 3.1 Poisson use model 

193 A natural way to motivate this specific model of space usage is to assume that individuals 

194 make a sequence of random resource selection decisions so that the outcomes n(x) (for all 

195 x) are marginally independent Poisson random variables: 



197 In this case, the number of visits to any particular cell is affected by the covariate z(x) 

198 but has a baseline rate (exp(ao)) related to the amount of movement occurring over some 

199 time interval. This is an equivalent model to the multinomial model given previously in the 

200 sense that, if we condition on the total sample size R = J2x^i^)y then the vector of use 

201 frequencies {n(x)} for individual with activity center s, has a multinomial distribution with 

202 probabilities 



ra(x)|s ~ Poisson(A(x|s)) 



196 



where 
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203 which is the same as Eq. [2] because cancels from the numerator and denominator of the 

204 multinomial cell probabilities and thus this parameter is not relevant to understanding space 

205 usage. Note that if use frequencies are summarized over i = 1,2,..., Nt^i individuals for 

206 each pixel, then a standard Poisson regression model for the resulting "quadrat counts" is 



207 reasonable. This corresponds to "Design 1" in Manly et al. (2002). 



208 3.2 Random Thinning 

209 Suppose our sampling is imperfect so that we only observe a smaller number of telemetry 

210 fixes than actual use frequency, ?7,(x). We express this "thinning" (or sampling) by assuming 

211 the observed number of uses is a binomial random variable based on a sample of size ?t,(x): 

m(x) ~ Bin(n(x), 0o)- 

212 Then, the marginal distribution of the new random variable m is also Poisson but with 

213 mean 

log{X{x\s)) = log{(f)o) + ao — ai(i(x|s)^ + a2z{'x). 

214 Thus, the space-usage model (RSF) for the thinned counts m is the same as the space-usage 

215 model for the original variables n. This is because if we remove n from the conditional model 

216 by summing over its possible values, then the vector of thinned use frequencies m (i.e., for 

217 all pixels) is also multinomial with cell probabilities 

I . A(x|s) 
■7r(x|s) — 



E.A(x|s) 
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218 and so the constants Oq and 0o cancel from both the numerator and denominator. Thus, 

219 the underlying RSF model applies to the true unobserved count frequencies n and also those 

220 produced by a random thinning or sampling process, m. 

221 In summary, if we conduct a telemetry study of i = 1, 2 . . . , Ntei individuals, the observed 

222 data are the nG x 1 vectors of use frequencies riij for each individual. We declare these data to 

223 be "resource-selection data" which are typical of the type used to estimate resource-selection 



224 functions (RSFs) (Manly et al. , 2002). In fact, the situation we have described here in which 



225 we obtain a random sample of use locations and a complete census of available locations is 



226 referred to as "Design 2" by (Manly et al. , 2002) 



227 3.3 Resource Selection in SCR Models 

228 The key to combing RSF data with SCR data is to work with this underlying resource 

229 utilization process and formulate SCR models in terms of that process. Imagine that we have 

230 a sampling device, such as a camera trap, in every pixel. If the device operates continually 

231 then it is no different from a telemetry instrument. If it operates intermittently or does not 

232 expose the entire area of each pixel then a reasonable model for this imperfect observation is 

233 the "thinned" binomial model given above, where 0o represents the sampling effectiveness of 

234 the device. For data that arise from SCR studies, the frequency of use for each pixel where 

235 a trap is located serves as an intermediate latent variable that we don't observe. From a 

236 design standpoint, the main difference between SCR studies and telemetry is that, for SCR 

237 data, we do not have sampling devices in all locations (pixels) in the landscape. Rather, the 

238 data are only recorded at a subsample of them, the trap locations, which we identify by the 

239 specific coordinates Xi, . . . , xj. 
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240 So we imagine that the hypothetical perfect data from a camera trapping study are the 

241 counts m(x) only at the specific trap locations x^, and for all individuals in the population 

242 i = 1,2, . . . , N where N > Ntei- We denote the individual- and trap-specific counts by rriij 

243 for individual with activity center Sj and trap location x^. In practice, many (perhaps most) 

244 of the rriij = m(xj,Sj) frequencies will be 0, corresponding to individuals not captured in 

245 certain traps. We then construct our SCR encounter probability model based on the view 

246 that these frequencies m^- are latent variables. In particular, under the SCR model with 

247 binary observations, we observe a random variable yij = 1 if the individual i visited the pixel 

248 containing trap j and was detected. We imagine that yij is related to the latent variable rriij 

249 being the event rriij > 0, as follows: 

Uij ~ Bern(pij) 

250 where 

Pij = PY{mij > 0) = 1 - exp(-A(xj|sj)) 

251 This is the complementary log-log link relating pij to log(Aij), setting Ajj = A(xj|si): 

cloglogipij) = log{\ij) 

252 where 

log(Ajj) = log(0o) + oo - aid{yLj,Sif + a2z{yij). 

253 and we collect the constants so that ckq = Zog'(0o) + Oo is the baseline encounter rate which 

254 includes the constant intensity of use by the individual and also the baseline rate of detection. 
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255 conditional on use. 



25a 4 The Joint RSF/SCR Likelihood 

257 To construct the likelihood for SCR data when we have auxiliary covariates on space 

258 usage or direct information on space usage from telemetry data, we regard the two samples 

259 (SCR and RSF) as independent of one another. In practice, this might not always be the 

260 case but (1) the telemetry data often come from a previous study; (2) Or, the individuals 

261 are not the same, or cannot be reconciled, even if telemetry study occurs simultaneously; 

262 (3) In cases where we can match some individuals between the two samples, regarding them 

263 as independent should only entail a minor loss of efficiency because we are disregarding 

264 more precise information on a small number of activity centers. Moreover, we believe, it is 

265 unlikely in practice to expect the two samples to be completely reconcilable and that the 

266 independence formulation is the most generally realistic. 

267 Regarding the two data sets as being independent, our approach here is to form the 

268 likelihood for each set of observations as a function of the same underlying parameters 

269 and then combine them. In particular, let >Cscr(Q;o, cki, CK2, -^; Yscr) be the likelihood for the 

270 SCR data in terms of the basic encounter probability parameters and the total (unknown) 

271 population size N, and let >Cj.s/(q;i, Q;2; iUrs/) be the likelihood for the RSF data based on 

272 telemetry which, because the sample size of such individuals is fixed, does not depend on N . 

273 Assuming independence of the two datasets, the joint likelihood is the product of these two 

274 pieces: 

scr 
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275 In what follows, we provide a formulation of each likelihood component. An R function for 

276 obtaining the MLEs of model parameters is given in Appendix 1. 

277 We adopt the notation /(■) to indicate the probability distribution of whatever observ- 

278 able quantity is in question, e.g., f{u) is the marginal distribution of u and f{u\v) is the 

279 conditional distribution of u given v, etc. We use (?(•) to represent the probability distribu- 

280 tion of latent variables. The observation model for the SCR data for individual i and trap 

281 j, from sampling over K encounter periods, is: 



(3) 



where 



Pij = p{d{yij, Si), z(xj); ol) = I - exp(-A 



and 



284 We emphasize that this is conditional on the latent variables Sj (which appear in the distances 

285 dij). For these latent variables we adopt the standard assumption of uniformity, Sj ~ Unif(iS) 



286 for each individual i = 1,2, . . . ,iV (Royle and Young, 2008) where S is the state-space of 

287 the random variable s. 

288 The joint distribution of the data for individual i, conditional on Sj, is the product of J 

289 binomial terms (i.e., the contributions from each of J traps): 



f{yi\si,a) = Y[Bm{K,pij{(x)). 

i=i 



15 



290 
291 
292 



The marginal likelihood (Borchers and Efford, 2008) is computed by removing Sj, by inte- 
gration, from the conditional-on-s likelihood and regarding the marginal distribution of the 
data as the likelihood. That is, we compute: 



f{yi\a)= / f{yi\Si,a)g{Si)dSi 



s 



293 where, under the uniformity assumption, we have g{s) = The joint likelihood for all 

294 individuals, is the product of N such terms: 



295 In practice, we don't know and so we can't just compute the SCR likelihood in this manner. 

296 Instead, we compute the contributions of the n observed individuals directly as given above, 

297 but then we have to compute the likelihood contribution for the "all 0" encounter history, 

298 i.e., that corresponding to unobserved individuals. The mechanics of computing that are the 

299 same as for an ordinary observed encounter history, requiring that we integrate a binomial 

300 probability of over the state-space S: 



301 We then have to deal with the issue that n itself is a random variable, and that leads to 

302 the combinatorial term in front of the likelihood which involves the total population size 



N 





303 



This produces the conditional-on-A^ or "binomial form" of the likelihood (Borchers and 
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Efford, 


2008 


Royle , 


2009 



n\{N -n)\ 



N-n 




. i=l 



305 For the RSF data from the sample of individuals with telemetry devices we adopt the 

306 same basic strategy of describing the conditional-on-s likelihood and then computing the 

307 marginal likelihood by averaging over possible values of s. We have nij, the nG x 1 vector 

308 of pixel counts for individual i, where these counts are derived from a telemetry study or 

309 similar. We index these elements as rriig for individual i and grid cell g, noting that our index 

310 j is reserved only for trap locations, which are a subset of the nG coordinates Xi, . . . jX^g. 

311 The conditional-on-Sj distribution of the telemetry data from individual i is, omitting the 

312 multinomial combinatorial term which does not depend on parameters, 

nG 

f{mi\si,a) oc Jj7r(x3|si)"''' 

9=1 



313 where 



7rl,Xg|SiJ 



314 The marginal distribution is 



/(m,;|Q;) = J f{uii\Si,cx)g{Si)dSi 



315 and therefore the likelihood for the RSF data is 



Ntel 



£^^/(Q;|mi, ms, . . . , ramei) = JJ /(mi|Q;) 



1=1 
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316 A key technical aspect of computing these hkehhoods is the evaluation of the 2-dimensional 

317 integral over the state-space S, which we approximate (Appendix 1) by a summation over a 

318 fine mesh of points. We note also that the binomial form of the likelihood here is expressed 

319 in terms of the parameter N, the population size for the landscape defined by S. Given S, 

320 density is computed as D{S) = N/aiea^S). In our simulation study below we report N as 



321 the two are equivalent summaries of the data once S is defined. Borchers and Efford (2008) 

322 develop a likelihood based on a further level of marginalization, in which N is removed from 

323 the likelihood by averaging over a Poisson prior for A^. 



324 5 Simulation Analysis 

325 We carried-out a simulation study using the landscape shown in Fig. [T| and based on 

326 populations of size A^ = 100 and A^ = 200 individuals with activity centers distributed 

327 uniformly over the landscape. This covariate was simulated by generating a field of spatially 

328 correlated noise to emulate a typical patchy habitat covariate relevant to habitat quality 

329 for a species. We subjected individuals to sampling over A' = 10 sampling periods, using 

330 a 7 X 7 array of trapping devices located on the the integer coordinates {u * 5,v * 5) for 

331 u,v = 1, 2, 3, 4, 5, 6, 7. The SCR encounter model was of the form 

cloglog(pij) = tto - + "22(xi) 

332 with ao = —2, a = 2 and 02 = 1- In the absence of the covariate z, this corresponds to a 

333 RSF that is bivariate normal with standard deviation 2. These settings yielded an average of 

334 about n = 61 individuals captured for the A^ = 100 case and about n = 123 for the A^ = 200 
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335 case. The latter case represents what we beheve is an extremely large sample size based on 

336 our own experience and thus it should serve to gauge the large sample bias of the likelihood 

337 estimator. 

338 In addition to simulating data from this capture-recapture study, we simulated 2, 4, 8, 12, 

339 1 6 telemetered individuals to assess the improvement in precision as sample size increases. 

340 For all cases we observed 20 telemetry fixes per individual, assuming individuals were using 

341 space according to a RSF model with the same parameters as those generating the SCR data. 

342 We simulated 500 data sets for each scenario and, for each data set, we fit 3 models: (i) the 

343 SCR only model, in which the telemetry data were not used; (ii) the integrated SCR/RSF 

344 model which combined all of the data for jointly estimating model parameters; and (iii) 

345 the RSF only model which just used the telemetry data alone (and therefore ao and N are 

346 not estimable parameters). The focus of the simulations was to address the following basic 

347 questions: (1) how much does the root mean-squared error (RMSE) of N improve as we add 

348 or increase the number of telemetered individuals? (2) How well does the SCR model do 

349 at estimating the parameter of the RSF with no telemetry data? (3) How much does the 

350 precision of the RSF parameter improve if we add SCR data to the telemetry data? 

351 Results for = 100, = 200 and Ntei = (2, 4, 8, 12, 16) are presented in Table |lj We 

352 note that the first row of each batch (labeled "SCR only" ) represent the same estimator and 

353 data configuration. These replicate runs of the SCR-only situation give us an idea of the 

354 inherent MC error in these simulations, which is roughly about 0.25 and 0.89 on the scale 

355 for the N = 100 and = 200 cases, respectively. The mean A^ for the SCR-only estimator 

356 across all 5 simulations for A^ = 100 was mean(iV) = 99.418, an empirical bias of 0.6%. For 

357 A^ = 200, the estimated A^ across all 5 simulations (5 levels of Ntei) was mean(A^) = 199.712, 
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358 an empirical bias of about 0.15%, within the MC error of the true value of = 200. The 

359 results suggests a very small bias of < 1% in the MLE of for both the SCR-only and 

360 combined SCR/RSF estimators. In practice, we expect a small amount of bias in MLEs as 

361 likelihood theory only guarantees asymptotic unbiasedness. 

362 In terms of RMSE for estimating A^, we see that (Table [T]), generally, there is about a 5% 

363 reduction in RMSE when we have at least 2 telemetered individuals. And, although there is a 

364 lot of MC error in the RMSE quantities, it might be as much as a 10% reduction as the sample 

365 size of captured individuals increases under the higher A^ = 200 setting. This incremental 

366 improvement in RMSE of N makes sense because, while the telemetry provides considerable 

367 information about the structural parameters of the model, it provides no information about 

368 mean p, i.e. ao, which comes only from the SCR data. Thus estimating A^ benefits only 

369 slightly from the addition of telemetry data. 

370 The MLE of the RSF parameter a2 exhibits negligible or no bias under both the SCR 

371 only and SCR/RSF estimators. It is well-estimated from SCR data alone and even better 

372 than RSF data alone (in terms of RMSE) until we have more than 200 or so telemetry 

373 observations. The biggest improvement from the use of telemetry data comes in estimating 

374 the parameter a. We see that a is effectively unbiased, and there is a very large improvement 

375 in RMSE of a, perhaps as much as 50-60% in some cases, when the telemetry data are used in 

376 the combined estimator (that really doesn't translate much into improvements in estimating 

377 A^ as we saw previously). Improvement due to adding telemetry data diminishes as the 

378 expected sample sizes increases, and so telemetry data does less to improve the precision of 

379 a and a2 for A^ = 200 than for A^ = 100. This is because the SCR data along are informative 

380 about both of those parameters. 
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381 The results as they concern hkehhood estimation of suggest that there is not a substan- 

382 tial benefit to having telemetry data. Estimators "SCR only" and "SCR/RSF" both appear 

383 approximately unbiased for = 100 and = 200, and for any sample size of telemetered 

384 individuals. The RMSE is only 5-10% improved with the addition of telemetry information. 

385 However, we find that there is substantial bias in N if we use the misspecified model that 

386 contains no resource selection component. That is if we leave the covariate z{x) out of the 

387 model and incorrectly fit a model with symmetric and spatially constant encounter model, 

388 we see about 20% bias in the estimates of in a limited simulation study that we carried- 

389 out (Tab. |2|. As such, accounting for resource selection is important, even though, when 

390 accounted for, telemetry data only improves the estimator incrementally. In addition, we 

391 find that the importance of telemetry data is relatively more important for smaller sample 

392 sizes. We carried-out one simulation study for the A^ = 100 case but with lower average 

393 encounter probabilities, setting ao = —3. This produces relatively smaller data sets with 

394 E[n\ = 37. The results are shown in Tab. |3} There are some important features evident 

395 from this table. First, as a result of the small samples, the MLE of A^ is biased for both SCR 

396 only and SCR/RSF estimators although less biased for the SCR/RSF estimator than for 

397 SCR only. The persistent bias in A^ for both models results from the information about ao 

398 coming only from SCR data, and that estimator itself is intrinsically biased in small samples. 

399 Conversely, the estimator of a2, the RSF parameter, appears unbiased for all 3 estimators 

400 (SCR only, SCR/RSF and RSF only), as does the estimator of cr. We see relatively larger 

401 improvements in RMSE (compared with Tab. [l]) of N, and those improvements increase 

402 substantially as Ntei increases. 
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403 6 Discussion 



How animals use space is a fundamental interest to ecologists, and important in the 
conservation and management of many species. Normally this is done by telemetry and 



models referred to as resource selection functions (Manly et al. , 2002). Conversely, spatial 



capture- recapture models have grown in popularity over the last several years (Efford, 2004 
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2012). These, and indeed, most, development and applications of SCR models have focused 
on density estimation, not understanding space usage. However, it is intuitive that space 
usage should affect encounter probability and thus it should be highly relevant to density 
estimation in SCR applications. Despite this, a description of the relationship between 
encounter probability and space usage has not been developed in the literature on spatial 
capture-recapture models. Essentially all published applications of SCR models to date have 
been based on simplistic encounter probability models that are symmetric and do not vary 



417 across space. One exception is Royle et al. (2012) who developed SCR models that use 



ecological distance metrics ("least-cost path") instead of normal Euclidean distance. Here 
we developed an SCR model in terms of a basic underlying model of space or resource use. 



that is consistent with existing views of resource selection functions (RSFs) (Manly et al. 



2002j). 

In developing the SCR model in terms of an underlying model of space usage, we achieve 
a number of enormously useful extensions of existing SCR and RSF methods: (1) We have 
shown how to integrate classical RSF data from telemetry with spatial capture-recapture 
data based on individual encounter histories obtained by classical arrays of encounter devices 
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426 or traps. This leads to an improvement in our ability to estimate density, and also an 

427 improvement in our ability to estimate parameters of the RSF function. Thus, the combined 

428 model is both an extension of standard SCR models and also and extension of standard 

429 RSF models. As many animal population studies have auxiliary telemetry information, the 

430 ability to incorporate such information into SCR studies has enormous applicability and 

431 immediate benefits in many studies. While adding RSF data to SCR data may increase 

432 precision of the MLE of N only incrementally, the effect can be more substantial in sparse 

433 data sets and, generally, RSF produces relatively huge gains in precision in the MLE of a. (2) 

434 We have shown that one can estimate RSF model parameters directly from SCR data alone. 

435 While further exploration of this point is necessary, it does establish clearly that SCR models 

436 are explicit models of space usage. Because capture-recapture studies are, arguably, more 

437 widespread than telemetry studies alone, this greatly broadens the utility and importance of 

438 data from those studies. (3) It is also now clear that one of the important parameters of SCR 

439 models, that controlling "home range radius" , can be directly estimated from telemetry data 

440 alone. The combined RSF+SCR model does yield large improvements in estimation of a. As 

441 a practical matter, this suggests we could estimate a entirely from data extrinsic to the SCR 

442 study which might provide great freedom in the design of SCR studies. For example, traps 

443 could be spaced far enough apart to generate relatively few (even no) spatial recaptures, 

444 but dramatically increase the coverage of the population, i.e., the observed sample size of 

445 captured individuals relative to N. (4) Finally, we found that an ordinary SCR model with 

446 symmetric encounter probability model produces extremely biased estimates of when the 

447 population of individuals does exhibit resource selection. As such, it is important to account 

448 for space usage when important covariates are known to influence space usage patterns. 
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449 Use of telemetry data in capture-recapture studies has been suggested previously. For ex- 



450 ample, White and Shenk (2001) and Ivan (2012) suggested using telemetry data to estimate 



the quantity "probability that an individual is exposed to sampling" but their estimator 
requires that individuals are sampled in proportion to this unknown quantity, which seems 
impossible to achieve in many studies. In addition, they do not directly integrate the teleme- 
try data with the capture-recapture model so that common parameters are jointly estimated. 



In fact, they don't acknowledge shared parameters of the two models. SoUmann et al. (2012) 



did recognize this, and used some telemetry data to estimate directly the parameter a from 
the bivariate normal SCR model in order to improve estimates of density. This was an 
important conceptual development in the sense that it recognized the relationship between 
SCR models and models of space usage, but their model did not include an explicit resource 
selection component, and they did not implement a joint estimation framework. 



461 „ We developed a formal analysis framework here based on marginal likelihood (Borchers 



and Efford, 2008). In principle, Bayesian analysis does not pose any unique challenges for 



this new class of models although we expect some loss of computational efficiency due to 
the increased number of times the components of the likelihood would need to be evaluated. 
We imagine that some problems would benefit from a Bayesian formulation, however. For 
example, using an open population model that allows for recruitment and survival over time 



467 (Gardner et al. , 2010a) is convenient to develop in the BUGS language and incorporating 



468 information on unmarked individuals has been done using Bayesian formulations of SCR 



469 models (Chandler and Royle, 2012 SoUmann et al. , 2012) but, so far, not likelihood methods. 



In our formulation of the joint likelihood for RSF and SCR data, we assumed the data 
from capture-recapture and telemetry studies were independent of one another. This implies 
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472 that whether or not an individual enters into one of the data sets has no effect on whether 

473 it enters into the other data set. We cannot foresee situations in which violation of this 

474 assumption should be problematic or invalidate the estimator under the independence as- 

475 sumption. In some cases it might so happen that some individuals appear in both the RSF 

476 and SCR data sets. In this case, ignoring that information should entail only an incremental 

477 decrease in precision because a slight bit of information about an individuals activity center 

478 is disregarded. Heuristically, an SCR observation (encounter in a trap) is like one additional 

479 telemetry observation, and so the misspecification (independence) regards the two pieces of 

480 information as having separate activity centers. Our model pretends that we don't know 

481 anything about the telemetered individuals in terms of their encounter history in traps. In 

482 principle it shouldn't be difficult to admit a formal reconciliation of individuals between the 

483 two lists. In that case, we just combine the two conditional likelihoods before we integrate 

484 s from the conditional likelihood. This would be almost trivial to do if all individuals were 

485 reconcilable (or none as in the case we have covered here) but, in general , we think you will 

486 always have an intermediate case - i.e., either none will be or at most a subset of teleme- 

487 tered individuals will be known. More likely you have variations of "well, that guy looks 

488 telemetered but we don't know which guy it is....hmmm" and that case, basically a type of 

489 marking uncertainty or misclassification, is clearly more difficult to deal with. 

490 We conclude that the key benefit of our combined SCR/RSF model is its ability inte- 

491 grate realistic patterns of space usage directly into SCR models and avoid extreme bias in 

492 estimating and, secondarily, we are able to obtain RSF information from SCR alone. 

493 Therefore, our new class of integrated SCR/RSF models allows investigators to model how 

494 the landscape and habitat influence movement and space usage of individuals around their 
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495 home range, using non-invasively collected capture-recapture data or capture-recapture data 

496 augmented with telemetry data. This should improve our ability to understand, and study, 

497 aspects of space usage and it might, ultimately, aid in addressing conservation-related prob- 

498 lems such as reserve or corridor design. And, it should greatly expand the relevance and 

499 utility of spatial capture-recapture beyond simply its use for density estimation. 
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572 ### before rvuming this code, put the functions at the end of this script 
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573 ### into your R workspace 

574 ### 

575 
576 

577 ## the following block of code makes up a covariate as a spatially correlated 

578 ## noise field, with an exponential spatial correlation function 

579 set.seed(1234) 

580 gr<-expand.grid(l :40, 1 :40) 

581 Dmat<-as. matrix (dist(gr)) 

582 V<-exp(-Dmat/5) 

583 z<-t (chol (V) )%*%rnorm(1600) 

584 spatial, plot (gr,z) 

585 
586 

587 ### 

588 ### Set some parameter values 

589 ### 

590 alphaO <- -2 

591 sigma<- 2 

592 beta<- 1 

593 Ntel<-4 # number of individuals with telemeters 

594 nsim<-100 

595 Nfixes<-20 # number of telemetry fixes per individual 

596 N<- 100 # population size 

597 
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598 

599 # simulate activity centers of all N individuals 

600 Sid<- sampled : 1600, N,replace=TRUE) 

601 # and coordinates 

602 S<-gr[Sid,] 

603 # now draw centers of telemetered individuals 

604 # have to draw telemetry guys interior or else make up more landscape — 

605 # can't have truncated telemetry obs 

606 

607 poss.teK- S[,l]>5 & S[,l]<35 & S[,2]>5 & S[,2]<35 

608 tel . guys<-sample (Sid [poss . tel] , Ntel) 

609 sid<-tel . guys 

610 stel<-gr [sid,] 

611 

612 # make up matrix to store RSF data 

613 n< -matrix (NA , nrow=Ntel , ncol=1600) 

614 

615 # for each telemetered guy simulate a niimber of fixes. 

616 # note that n = for most of the landscape 

617 par(mfrow=c(3,3) ) 

618 lammat<-matrix (NA , nrow=Nt el , ncol=1600) 

619 for(i in l:Ntel){ 

620 d<- Dmat [sid[i] ,] 

621 lam<- expd - (l/(2*sigma*sigma) ) *d*d + beta* z) 

622 n [i , ] <-rmult inom ( 1 , Nf ixes , lam/ sum (lam) ) 
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623 par(mar=c(3,3,3,6)) 

624 lammat [i,] <-lam 

625 img<- matrix (lam, nrow=40,ncol=40,byrow=FALSE) 

626 imaged : 40, 1 : 40, rot (img) , col=terrain. colors (10) ) 

627 3" 

628 

629 ## now lets simulate some SCR data on a bunch of guys: 

630 

631 # make a trap array 

632 X<- cbind( sort(rep( seq(5,35,5) ,7)) , rep( seq(5,35,5) ,7)) 

633 ntraps<-nrow(X) 

634 raster .point<-rep(NA,nrow(X) ) 

635 for(j in l:nrow(X)){ # which piont in the raster is the trap? must be raster points 

636 raster.point[j]<- (1:1600) [ (X[j , l]==gr [, 1] ) & (X[j,2] == gr[,2])] 

637 y 

638 points (X,pch=20,cex=2) 

639 

640 D<- e2dist(S,X) ## N x ntraps 

641 Zmat<- matrix (z [raster .point] ,nrow=N,ncol=ntraps,byrow=TRUE) # note make dims the same 

642 loglam<- alphaO -(l/(2*sigma*sigma) )*D*D + beta*Zmat 

643 p<- l-exp(-exp(loglam) ) 

644 

645 ## Now simulate SCR data 

646 

647 K<- 10 
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648 y<-matr ix (NA , nrow=N , ncol=ntraps) 

649 for(i in 1:N){ 

650 y[i,]<- rbinom(ntraps,K,p[i,] ) 

651 3" 

652 

653 cap<-apply(y , 1 , suin)>0 

654 

655 y<-y[cap,] 

656 gr<-as .matrix (gr) 

657 sbar<- (nyo*yogr) /as. vector (n°/o*°/orep(l,nrow(gr))) 

658 

659 # Basic SCR model with RSF covariate at trap locations. 

660 tmpl<-nlm(intlik3rsf . v2,c(-3,log(3) ,1,0) ,y=y ,K=K,X=X,ztrap=z [raster .point] ,G=gr) 

661 

662 # use telemetry data and activity centers for those are marginalized out of the likelihood 

663 tmp2<-nlm ( intlikSrsf . v2 , c (-3 , log (3) , 1 , 0) , y=y , K=K , X=X , ztrap=z [raster . point] , G=gr , ntel=n , zall=as 

664 

665 # use mean "s" instead of estimating it 

666 tmp3<-nlm ( intlik3rsf . v2 , c (-3 , log (3) , 1 , 0) , y=y , K=K , X=X , ztrap=z [raster . point] , G=gr , ntel=n , zall=as 

667 

668 # no SCR data, s is random. Here there are 2 extra parameters that are not estimated: start [1] 

669 tmp4<-nlm(intlik3rsf . v2,c(-3,log(3) ,1,0) ,y=NULL,K=K,X=X,ztrap=z [raster. point] ,G=gr ,ntel=n,zall= 

670 

671 # Fits SCR model with isotropic Gaussian encounter model 

672 tmp5<- nlm(intlik3rsf . v2,c(-3,log(3) ,1,0) ,y=y,K=K,X=X,ztrap=rep(0,ntraps) ,G=gr) 
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673 
674 

675 ### 

676 ### put all the functions below this line into your R workspace 

677 ### 

678 

679 spatial .plot<- 

680 function(x,y){ 

681 nc<-as. numeric (cut (y, 20)) 

682 plot(x,pch=" ") 

683 points (x,pch=20,col=topo. colors (20) [nc] ,cex=2) 

684 ###image . scale (y , col=topo . colors (20) ) 

685 3" 

686 

687 ### This is the likelihood function 

688 ### It computes several versions of the likelihood depending on the arguments specified 

689 ### see the 5 examples above 

690 

691 int likSr sf . v2 <-f unction (start=NULL , y=y , K=NULL , X=traplocs , ztrap , G , ntel=NULL , zall=NULL , stel=NUL] 

692 # start = vector of length 5 = starting values 

693 # y = nind x ntraps encounter matrix 

694 # K = how many samples? 

695 # X = trap locations 

696 # ztrap = covariate value at trap locations 

697 # zall = all covariate values for all nG pixels 
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698 # ntel = nguys x nG matrix of telemetry fixes in each nG pixels 

699 # stel = home range center of telemetered individuals, IF you wish to estimate it. Not 

700 

701 nG<-nrow(G) 

702 D<- e2dist(X,G) 

703 

704 alphaO<-start [1] 

705 sigma<- exp(start [2] ) 

706 alpha2<- start [3] 

707 nO<- expCstart [4] ) 

708 aO<- 1 

709 

710 if ( ! is .null(y) ){ 

711 loglam<- alphaO -(l/(2*sigma*sigma))*D*D + alpha2*ztrap # ztrap recycled over nG 

712 probcap<- l-exp(-exp(loglam) ) 

713 #probcap<- (exp(thetaO)/(l+exp(thetaO)))*exp(-thetal*D*D) 

714 Pm<-matr ix (NA , nrow=nrow (probcap) , ncol=ncol (probcap) ) 

715 ymat<-y 

716 ymat<-rbind(y,rep(0,ncol(y))) 

717 lik.marg<-rep(NA,nrow(ymat) ) 

718 for(i in l:iirow(ymat)){ 

719 Pm[l: length (Pm)]<- (dbinom(rep(ymat [i ,] ,nG) ,rep(K,nG) .probcap [1 : length (Pm)] ,log=TRUE)) 

720 lik.cond<- exp(colSums(Pm)) 

721 lik.marg[i] <- sum( lik.cond*(l/nG) ) 

722 3" 
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723 nv<-c(rep(l,length(lik.inarg)-l) ,iiO) 

724 partK- lgainina(nrow(y)+nO+l) - IgammaCnO+l) 

725 part2<- sum(nv*log(lik.inarg)) 

726 out<- -l*(partl+ part2) 

727 3" 

728 else-[ 

729 OUt<-0 

730 3" 

731 

732 if (!is.null(ntel) & ! is .null(stel) ){ 

733 

734 # this is a tough calculation here 

735 D2<- e2dist(stel,G)"2 

736 # lam is now nG x nG! 

737 lam<- t(exp(aO - (l/(2*sigma*sigma))*t(D2)+ alpha2*zall) ) # recycle zall over all ntel guys 

738 denom<-rowSums(lam) 

739 probs<- lam/denom # each column is the probs for a guy at column [j] 

740 

741 tel.loglik<- -l*sum( ntel*log (probs) ) 

742 

743 out<- out + tel.loglik 

744 3" 

745 

746 if (! is .null (ntel) & is .null(stel) ){ 

747 
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748 # this is a tough calculation here 

749 D2<- e2dist(G,G)"2 

750 # lam is now nG x nG! 

751 lam<- t(exp(aO - (l/(2*sigma*sigma))*t(D2)+ alpha2*zall) ) # recycle zall over all ntel guys 

752 denom<-rowSums(lam) 

753 probs<- t (lam/denom) # each column is the probs for a guy at column [j] 

754 marg<- as.vector(rowSums(exp(nteiyo*%log(probs))/nG )) 

755 

756 tel.loglik<- -l*sum(log(marg) ) 

757 

758 out<- out + tel.loglik 

759 y 

760 

761 out 

762 y 

763 
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Table 2: Expected value of N and a for truth N — 200 and a — 2 under a model of resource 
selection with a single covariate, when the encounter probability model is misspecified by 
a symmetric and constant model assuming no resource selection; column "bias" is percent 
bias. 







E[N] 


bias 


RMSE 


E 


a 


RMSE 


n= 


=2 


161.48 


-19.2 


39.98 


1.84 


0.180 


n= 


=4 


161.32 


-19.3 


40.00 


1. 


83 


0.191 


n= 


=8 


161.46 


-19.3 


40.06 


1. 


84 


0.184 


n= 


12 


162.40 


-18.8 


38.95 


1. 


84 


0.185 


n= 


16 


160.93 


-19.5 


40.44 


1.84 


0.190 



Table 3: Mean and RMSE of the 
for the = 100 and "low p" case, 
the SCR likelihood only, the joint 
the latter, the parameter N is not 



sampling distribution of the MLE for model parameters 
For each of 500 simulated data sets, a model was fit using 
SCR/RSF likelihood, and the RSF hkelihood only. For 
statistically identifiable. 



Estimator E[N] RMSE E[a2] RMSE E 



a 



RMSE 



Ntel = 2 






















SCR only 
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22.88 
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00 





19 
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02 





261 


SCR/RSF 
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20.98 
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00 





17 


2 


00 





136 


RSF only 
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02 





30 


1 


99 





163 


Ntel = 4 






















SCR only 
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26.52 
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01 
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01 





258 


SCR/RSF 


103.55 


22.92 
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01 





14 
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00 
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RSF only 
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01 





21 


1 


99 





114 


Ntel = 8 






















SCR only 


107.41 


45.05 
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19 
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SCR/RSF 
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22.13 
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00 
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00 





076 


RSF only 
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1 


99 





081 


Ntel = 12 






















SCR only 
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RSF only 
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RSF only 
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00 
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FIGURE CAPTIONS 



765 Figure 1: A typical habitat covariate reflecting habitat quality or hypothetical utility of 

766 the landscape to a species under study. Home range centers for 8 individuals are shown with 

767 black dots. 

768 Figure 2: Space usage patterns of 8 individuals under a space usage model that contains a 

769 single covariate (shown in Fig. [T]). Plotted value is the multinomial probability vTjj for pixel 

770 j under the model in Eq. [2j 
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Figure 1: A typical habitat covariate reflecting habitat quality or hypothetical utility of the 
landscape to a species under study. Home range centers for 8 individuals are shown with 
black dots. 
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Figure 2: Space usage patterns of 8 individuals under a space usage model that contains a 
single covariate (shown in Fig. [T]). Plotted value is the multinomial probability for pixel 
j under the model in Eq. [2] 
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